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We demonstrate that dynamical symmetry plays a crucial role in determining the structure of the 
eigenspectra of spinor Bose-Einstein condensates (BECs). In particular, the eigenspectra of spin- 
1 and spin-2 BECs in the single-mode approximation are shown to be completely determined by 
dynamical symmetries, where a spin-2 BEC corresponds to the (7(5) limit of the interacting boson 
model in nuclear physics. The eigenspectrum of a spin-3 BEC is solved analytically for a specific 
class of coupling constants, while it is shown that dynamical symmetry alone is not sufficient to 
determine the spectrum for arbitrary coupling constants. We also study the low-lying eigenspectra 
of spin-1 and spin-2 BECs in the absence of external magnetic field, and find, in particular, that 
the quasi-degenerate spectra emerge for antiferromagnetic and cyclic phases. This implies that 
these systems are highly susceptible to external perturbations and may undergo symmetry-breaking 
transitions to other states upon increasing the system's size. 

PACS numbers: 03.75.Hh, 03.75.Mn, 05.30.Jp 



Eigenvalue problems that can be solved analytically give profound insight into quantum many-body systems and 
are also important to elucidate the nature of the problems. In many of such cases, the solvability originates from the 
symmetry of a system. Here, the symmetry includes not only the space-time symmetry but also dynamical symmetry 
which arises from the special properties of the forces. If a system has a dynamical symmetry, the eigenvalue problem 
can be solved algebraically and we can find the exact spectrum which reveals the structure of the system. For instance, 
collective spectra of a number of atomic nuclei can be understood from dynamical symmetries [ij . 



The U{n) symmetry in the n-dimensional harmonic-oscillator problem and the 0(4) symmetry in the three- 
dimensional Coulomb problem are two well-known examples of dynamical symmetries [2|. Dynamical symmetries 
have played crucial roles in such diverse fields of physics as elementary particle physics p\, nuclear physics and 
molecular physics [4j]. In the context of Bose-Einstein condensates (BECs) of dilute atomic vapor, the S'[/(l, 1) sym- 
metry in a harmonically-trapped two-dimensional system was discovered by Pitaevskii and Rosch Q . In the present 
paper, we analyze the dynamical symmetry of spinor BECs. 

Mean-field theories of spinor BECs have previously been discussed for spin-1, 2, and 3 cases in, e.g., Refs. @, 
d, H, and dH [H, [11] , respectively. The ground-state phases of spin-1 BECs comprise ferromagnetic and polar 
phases. In addition to these two phases, the ground-state phases in a spin-2 BEC involves a cyclic phase which can 
hold 1/3-quantum vortices [13, [HI. A rich variety of phases have been predicted in spin-3 BECs 11, 12, 13, 20, 2lj; 
however, their physical properties have yet to be fully investigated. 

The exact eigenspectra and eigenstates have been obtained for spin-1 and spin-2 BECs [1, [l^, [H, [IE [III- For 
the case of antiferromagnetic coupling, it is found that the ground state of a spin-1 BEC is a spin-singlet pair-boson 
condensate which is frag mented in the sense that more than one eigenvalue of the single-particle density operator is 
of the order of I, [M [H . As a consequence, the number of atoms in each magnetic sublevel fluctuates violently 
between 1 and A^ as the magnetization of a system increases [3, [l^. For the case of a spin-2 BEC, on the other 
hand, the cyclic and antiferromagnetic phases are predicted to exhibit a spin-singlet trio-boson condensate, and a 
Meissner-like effect in the magnetic response, respectively [1, [l3|- 

In this paper, we develop a systematic method of representing the eigenspectra and eigenstates of a spinor BEC 
by exploiting the dynamical symmetry of the system, and show that the eigenspectra and eigenstates of spin-1 and 
spin-2 BECs are determined for the entire range of coupling constants and those of spin-3 BEC are determined for 
a specific class of coupling constants. In particular, a spin-2 BEC corresponds to the U{h) limit of the interacting 
boson model, which describes the collective properties of atomic nuclei [l|, and the U{b) classification can be utilized 
for representing excitation spectra of a spin-2 BEC. We examine the low- lying eigenspectra and eigenstates of spin-1 
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and spin-2 BECs in the absence of external magnetic field and find, in particular, that the quasi-degenerate spectra 
emerge above the ground states of antiferromagnetic and cyclic phases. This fact suggests that these ground states 
are highly susceptible to external perturbations and may undergo symmetry-breaking transitions to other states as 
the size of the system increases. 

This paper is organized as follows. Section It presents the Hamiltonian of a spin-/ BEC in terms of spherical tensor 
operators. Section IH briefiy reviews elements of group theory which are relevant to later discussions. Sections IV 
, V, and VI use the formalism described in Sec. HI to analyze spin-1, -2, and -3 BECs, respectively. Section VII 
investigates the low- lying eigenspectra and eigenstates of spin-1 and spin-2 BECs in the absence of external magnetic 
field. Section VIII summarizes the main results of the present paper. 



II. FORMULATION OF THE PROBLEM 



We consider a system comprised of identical bosons with spin / and mass M interacting via s-wave contact inter- 
action. The spin-independent part of the interaction is by far the largest and determines the density distribution of 
the particles. We shall therefore assume that all particles share a single spatial mode (single-mode approximation), 
and focus on the many-body spectrum of the spin state. 

The spin-dependent part of the interaction Hamiltonian can, in general, be written as 



nxi-X2) = 5(xi-X2) ^ gpVF, (1) 

F=0,2,...,2/ 

where gp = 'inh^ap/M is the coupling constant with ap being the s-wave scattering length in the total spin F channel, 
and Vp is the projection operator that projects the wave function of a pair of interacting atoms into the total spin F 
channel. The corresponding second-quantized Hamiltonian is given by 



F=0.2,...,2f M=-F 



Jd^< fmfn\FM >< FM\fm' fn > ^1^1^^'^^', (2) 



where < fmfn\FM > is the Clebsch-Gordan coefficient, and >I'm(x) (m = — /, — ...,/) represents the annihilation 
operator of a boson at position x with magnetic quantum number m. In Eq. ([2]), repeated indices (to, to ,n,n) are 
assumed to be summed. In the presence of a magnetic field, the Hamiltonian also includes the Zeeman term 

dlC^l{F'-)„,n^n, (3) 

where i^^ is the z component of the spin operator, and p is the product of the gyromagnetic ratio and the external 
magnetic field which is assumed to be applied in the z direction. In the following discussions, we shall write the sum 
of V and the Zeeman term as H and simply call it the Hamiltonian. In the single-mode approximation the field 
operator takes the form of 

^'„j ~ 0(x)a„,, (4) 

where (/'(x) is the spatial mode into which all bosons are assumed to condense, and a™ is the corresponding annihilation 
operator with to being the magnetic quantum number. 

The dynamical symmetry of the Hamiltonian can be analyzed in a transparent manner in terms of spherical tensor 
operators. The spherical tensor operators of rank /, T^, transform as irreducible tensors under rotation and satisfy 
the following commutation relations: 

[F,,fl]=mfl, (5) 

, Tl] = ^fif+l)-m{m±l)fi^,, (6) 
where = ^ TOaJ„a,„, and F± V f{f + 1) - iTi{m ± l)a|„_,_ia,„. 

m m 

We introduce a hierarchy of spherical tensor operators and tensor products as follows. We first note that operators 

hi, a™ = (-l)^'-"a_™, (7) 
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satisfy conditions ([5]) and ([6]) and are therefore spherical tensor operators of rank /. In terms of these operators, we 
can introduce spherical tensor products of rank I as 

We note that [a^ x (a* = ^ ' ^ 1, ■ ■ ■ satisfy conditions ([5|) and ([6]) with / = Z and are therefore the spherical 
tensors of rank I. Since the Hamiltonian is a scalar, it should be expressed in terms of scalar quantities which can be 
constructed from Eq. ([5]), as 



la'xa'l X \a X a\ \ = 2^ 1 1" x a' |„ x \a x a 



Jo 



= ^ <l^l^i\00>< fmfm\lfi><lfi\fnfn > al^al^rano,^'- (9) 

m^m ,n.n 

In fact, it can be shown that the Hamiltonian is expressed in terms of the spherical tensor products as 

2n 



^ = X7 ^ 9F J2 < frnfn\FM ><FM\fm fn > a)„a)^a^'k^' -p{F^)rnna)^^ 



F=0.,2,...,2f M=-F 



(10) 



F=0,2,...,2/ 



where $7 = (/ dxl^l"*) ^ is an effective volume. The last expression in Eq. (fTO|) can be obtained as follows. We use 

< FM\fmfn >=< F, -M\f, -to, /, -n > and < FMFM |00 >= ^ ^ Sj^j together with Eq. jT]) to show 

\'2F + 1 

that 

< FMlfmfn > a,^,a„, = (-1)™'+"' < FM\f,~m'j,~n > h„/a^, 

^ < F,^M\fm'fn > 

= V2F + 1 < F, A/, F, -Af|00 >< F, ^M\fm' fn > a^>a„> (11) 
Substituting this into the first lines of Eq. pO)) and comparing the result with Eq. we obtain the desired 

expression. In rewriting the Zeeman term, we have used < f, to, f, — mllO >= m(— 1)-'^^'". / — -, -. 

s J, \ y ) y /(2/ + l)(/ + l) 

The eigenvalue problem of Hamiltonian (|10|) thus reduces to expressing it in terms of invariant quantities that 
represent the underlying dynamical symmetry of the system. The invariant quantities can be found by using group 
theory which is described in the following section. 



III. ELEMENTS OF GROUP THEORY 



To make this paper self-contained, we briefly review some elements of group theory that are relevant to our analysis 
of spinor BECs [2]. The one-particle properties of a spin-/ BEC are described by the generators of the unitary group 
[7(2/ + 1), aj^a„, which obey the commutation relations 

[aL«ni a^ai^] = (^n^aLai^ - (Jmiyoj^an- (12) 

In addition, Bose symmetry requires that the bases for a system of A^-identical bosons constitute a totally symmetric 
irreducible representation. 



A. Racah form 



The Racah form provides an alternative representation for the generators of the unitary group U{2f + 1) and is 
suitable for describing the interaction Hamiltonian with rotational invariance. The Racah form can be constructed 
by the replacement of the generator a3„d„ of the U{2f + 1) group with its rotationally covariant form [d^ x d]^ 

= 0, • • • , 2/). The commutation relations of the latter are given by 
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[a^ X a]l^,[a^ X a]l^, = ^ ^ {21 + l){2r + I) < Iml m \l" m" > 

I ,ni 



I I I 
f f f 



(-1)' -(-l)'+'l[atxaL„, (13) 



where l^^j, ^ J J Wigner 6-j symbol whose definition and fundamental properties are given in Appendix. From 

the commutation relations (12), we find that [cl' x al^ with I = 1, 3, • • • , 2/ — 1 form a subgroup of the U{2f + 1) 



group which is referred to as the SO{2f + 1) group 

The number operator and the angular momentum operator are expressed in terms of spherical tensor products as 



N = ^2/ + l[at X a]° (14) 

and 



respectively, where Fi = = F.j., + iFy, F_i = F^ = F^ ~ iFy, and Fq — F^. 

The eigenvalues of the spinor Hamiltonian are related to the invariants called Casimir operators that commute with 
all the generators of the group. The Casimir operators relevant to spinor BECs are those of the C/(2/+ 1), 50(2/ + 1), 
and 50(2) groups. 

We first consider the linear Casimir operators of the U{2f + 1) and 50(2) groups. The Casimir operator that 
commutes with all the generators of the J7(2/ + 1) group is the number operator: 



C'i(C/(2/ + l)) = V27TT[a^xa]" 

= N. (16) 

The 50(2) group represents the rotation about the z-axis, and hereby has only one generator, with its linear Casimir 
operator given by 



C.(S0,2)). /ffl±iffi£±i)|a..ai; 

= (17) 

We next consider the quadratic Casimir operators of the U{2f + 1) and 50(2/ + 1) groups. The quadratic Casimir 
operators are constructed by taking scalar products of the two generators. The scalar product of tensor operators T' 
and t/' of rank I is defined by 

f -11^ = (-l)'\/2r+T[f ' X U^]". (18) 
The quadratic Casimir operator of the U(2f + 1) group is given by 

2/ 

02(1/(2/ + 1)) = ^[flt X a]' • [flt X 
1=0 

= N{N + 2f). (19) 
The quadratic Casimir operator of the 50(2/ + 1) group is given for / = 1 by 

02(50(3)) ^ /(/ + 1K2/ + I) [at X 5]i • [at x (20) 



[at X al' • [flt X al' 



and for / ^ 1 by 

02(50(2/ + 1)) = 2 

Z=l,3,...,2/-1 

= N{N + 2/ - 1) - (2/ + l)[[at X 0^]° x [a x 5]°]°, (21) 
where a numerical factor of 2 in Eq. (|2ip is introduced for the sake of convenience. 
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B. Branching problem and dynamical symmetry 

To uniquely characterize the basis set of the problem, we need to find a complete set of quantum numbers. This is 
done by introducing a complete chain of subgroups: 

GdG' dG" D ■■■ . (22) 

To find the complete chain of subgroups is the aim of the branching problem, which can be solved according to the 
well-established rules 0. Since the interaction Hamiltonian of spinor BECs has the SO{3) symmetry, we consider the 
following chain: 

U{2f + 1)d SO{2f + 1) D • • • D 5-0(3) D S0{2). (23) 

We also note that, to characterize the states uniquely, we must find 2/+ 1 quantum numbers in the totally symmetric 
representations of the U{2f + 1) group. 

When the Hamiltonian is expressed in terms of the Casimir operators of a chain of groups only, the system is said 
to possess dynamical symmetry, and the problem can be solved algebraically 0. In the following sections we will 
discuss the dynamical symmetry of spinor BECs. 



IV. SPIN-1 BEC 



We first apply the method presented in Sec. HI to the eigenvalue problem of a spin-1 BEC. Throughout this paper 
we consider only the s-wave scattering. The Bose symmetry then requires that the total spin of two colliding bosons 
be either or 2. The corresponding scalar quantities that appear in the interaction Hamiltonian can be obtained from 
Eqs. dm), (EII), and (jXI)) as 



[a^ X a^]° X [5 X d.]y = \n{N + 1) - \f^, 



(24) 



fa^ X atl2 X [a X dff 



1 



3\/5 



{2N^ -AN + F'^). 



Substituting Eqs. and into Eq. UHl), we obtain 

^ = ^ V2F+l5F[[a^ X ciY X [a X d.]^ - V2p[a^ x Sl] 



F=0,2 



^ [coN{N-l) + c,{F^-2N) 
coN{N + 2) + C2N + ciF^ 



1 

2n 



where the coupling constants Ci (i = 0, • • • ,2) are given by 



12 11 18 

Co = gffo + g52, ci = --.go + g.92, C2 = --go - -^92- 



(25) 



(26) 



(27) 



The last equality in Eq. ([26|) shows that the Hamiltonian is expressed as a sum of the Casimir operators of the J7(3), 
50(3), and 50(2) groups. This means that the Hamiltonian comprises the Casimir operators of the chain of the 
groups 



C/(3) D 50(3) D 50(2). 



(28) 



Since the Hamiltonian is expressed in terms of the Casimir operators of the chain of the groups only, the eigenvalue 
problem can be solved algebraically. In the spin-1 case, we must specify three quantum numbers to characterize the 
state uniquely, as mentioned in Sec. HI. In the present case, we can choose iV, F, and Fz as the desired quantum 
numbers which arise from the C/(3), 50(3), and 50(2) groups, respectively; therefore, we obtain the eigenstate as 



\N,F,Fz > 



(29) 
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The corresponding eigenvalue of the Hamiltonian is 

1 



2^ [coNiN - 1) + cFiF + 1)] - pF, 
2^[coiV(7V + 2) + C2N + ciF{F + 1)] - pF,, 



(30) 



where we use S"?/ (2) algebras in order to obtain the eigenvalue of F^ and the allowed values of F are F = N, N — 
2,N — 4, 1 or [l6j. The eigenvalue problem of a spin-1 BEC is thus completely solved by the dynamical symmetry 
of the system. It is worthwhile to note that while all the Casimir operators that appear in the Hamiltonian (j26p 
commute with each other, the Hamiltonian itself does not commute with the generators of U{3) group. This holds 
true for higher-spin cases to be described later. The present situation is thus quite different from the familiar examples 
of dynamical symmetries such as the n-dimensional harmonic-oscillator problem and the three dimensional Coulomb 
problem; in either of the latter two cases, the corresponding Hamiltonian commutes with the generators of the 
underlying group {U(n) or 0(4)). 

V. SPIN-2 BEC 

We next consider the case of a spin-2 BEC. The Bose symmetry requires that the total spin of two colliding bosons 
be 0, 2, and 4. The corresponding scalar quantities in the interaction Hamiltonian take the following forms: 



[[a^ X at]0 X [a X 5]"]" = i (iV^ + 37V - ($2(50(5))) , (31) 

[[at X at]2 X [a X = -1- f2N^ - 4N - F^ + 2C2{SO{5))) , (32) 

7v5 ^ ^ 

[[at X at]4 X [a X 'a]r - ^ (^^' ^^N+ ^P' - ^(72(50(5))) . (33) 
Substituting these relations into Eq. (|10p . we rewrite the Hamiltonian as 



F=0,2,4 

1 

20 
1 

m 



doN{N - 1) + di(F2 _ Qj^^ ^S+S- 

5 



-pF, 

dsNiN + 4) + d^N + d5C2{SO{5)) + diF'^] - pF, (34) 



where S+ — x at]", S- = S^_^, and the coupling constants di (z = 0, • • • ,5) are given by 

4 3,11, 10 3 

do = -g2 + yff4, di = --52 + ^34, d2= go - y52 + yff4, 

1 2 18 , 1 12 108 
dz = 7,90 + 3,92 + iT^gA, di = --go - —g2 - ^54, 

5 7 35 5 7 35 

12 3 

^5 = -g5o + y.92 - ■^54- (35) 

The last equality in Eq. ([M]) shows that the Hamiltonian comprises the Casimir operator of the J7(5), iS'0(5), 50(3), 
and 50(2) groups. The same situation arises in the C/(5) limit of the interacting boson model in nuclear physics [l|. 
As in the latter case, the Hamiltonian is expressed in terms of the Casimir operators of the chain of the groups 

C/(5) D 50(5) D 50(3) D 50(2), (36) 

and we can therefore find exact solutions, solely because of dynamical symmetry. In the spin-2 case, we need five 
quantum numbers to uniquely characterize the state as mentioned below Eq. 
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As in the spin-1 case, the quantum numbers N, F, and F^ arise from the (7(5), SO{3), and SO{2) groups, respectively. 
To find the quantum number of the 50(5) group, we consider the following commutation relations: 



[^„^±] = ±^±, [S+,S-] = -2S,, (37) 

ise operators i 

operators is expressed as 



where Sz = 2^ ^ 1' "^^^^^ operators satisfy SU{1, 1) commutation relations and the Casimir operator of these 



= -S+S- -Sz + si (38) 

The eigenvalues of S'^ and Sz are given by 

S^\S,Sz >^ S{S~l)\S,Sz >, Sz\S,Sz >=Sz\S,Sz >, 

(39) 

1 5 

where ^ ^ 4' follows from the relation 



S±\S, Sz >= y/S^z±Sz- S{S-1)\S, Sz > (40) 

that the eigenvalue of 5"^ for a given 5 takes on values 

Sz^ S,S+l,S + 2,--- . (41) 

We introduce a quantum number r, which denotes the number of particles that do not form spin-singlet pairs. Thus, 
S is expressed in terms of r by 

S^lr + l (42) 

and N takes on values 

iV = T,T + 2,T + 4,--- . (43) 

Since we can substitute and r for S and Sz, the eigenstates can be described by \S,Sz >= \N,t >. We thus find 
that the Casimir operator of the 5*0(5) group is expressed as 

02(50(5)) = iV(iV + 3) - 45+5_ 

= 45^- I (44) 

The eigenvalue of 0*2(50(5)) is then 

5 



02(50(5))|5,5, > = 



45(5-1)-^ 



S,Sz > 

= t{t + 3)\N,t> . (45) 

This means that the quantum number of the 50(5) group is t (t = iV, — 2, — 4, • • • , 1 or 0). We have thus 
determined four quantum numbers N,t,F, and Fz- However, yet another quantum number is needed to uniquely 
characterize the states. This additional quantum number is not directly related to the Casimir operators and is called 
the missing label Q . We can choose this quantum number as the number of spin-singlet trios of bosons [ld| . Thus 
the eigenstate of spin-2 BECs is given by 

|Af,T,n3o,F,F, >, (46) 

where 77.30 is the number of spin-singlet trios, F = A, A -I- 1, 2A — 2, 2A with r — 3n3o 4- A. (Note that 2A ~ 1 is 
mis sing .) The complete basis set and the relationship between the quantum numbers were pointed out in Ref. 
[23 , |23| . The corresponding exact eigenvalue of the Hamiltonian is given by 



^ 2n 



doN{N -1) + —{N^ t){N + r + 3) + di{F{F + 1) - 6N} 
5 



pFz 



= 2^ [d3N{N + 4) + diN + d5T{T + 3) + diF{F + 1)] - pFz. (47) 

The eigenvalue problem of a spin-2 BEC has thus been solved completely by dynamical symmetry alone. We note 
that the eigenspectrum is degenerate with respect to the missing label 77.30. 
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VI. SPIN-3 BEC 

We finally investigate the case of a spin-3 BEC. This is the case in which dynamical symmetry alone cannot, in 
general, solve the eigenvalue problem. We will explain the reason for this and determine a specific class of coupling 
constants for which the problem can be solved by dynamical symmetry alone. We begin by noting the following 
relations: 



55 3iril3iri53 
3 3 3(13 3 3(13 3 3 



0. (48) 



This means that [a^ x 5,]^,^ ^ and [a^ x 5]^^ ^ form a subgroup of the S0{7) group, which is referred to as the exceptional 
group G2 2]. The quadratic Casimir operator of the G2 group is given by 

C2{G2) = 2 ^ [at X 5]' • [at x d.]K (49) 
1=1,5 

Thus we can consider the following chain of groups: 

U{7) D SO{7) D G2 D SO{3) D S0{2). (50) 

The interaction Hamiltonian of spin-3 BECs has four independent coupling constants corresponding to the number 
of channels of binary collisions with total spin 0, 2, 4, and 6. Since the number of the quadratic Casimir operators in 
the chain of the groups is also four, one might expect that the Hamiltonian can be expressed in terms of the Casimir 
operators alone. Unfortunately, this is not the case. To show this, we note that the following identity can be derived 
from Eq. (|A.6P : 

[at X a]i • [at x a]^ - 2[at x a]^ • [at x a]^ + [at x a]^ • [fit x a]^ = 0. (51) 
From Eqs. (gS]), and ([HI]), we obtain 

C2{G2) - ^(72(50(7)). (52) 

Thus the number of independent quadratic Casimir operators is three rather than four, and dynamical symmetry is 
not sufhcient to completely determine the exact eigenvalues. If, on the other hand, we restrict ourselves to a specific 
class of the coupling constants, dynamical symmetry determines the eigenspectrum. In fact, using Eqs. (|2ip . 
and (jA.6[) . we obtain the following relations: 

[[at X at]" X [a X 5]"]° = i (iV^ + 57V - (72(50(7))) , (53) 
[[at X at]2 X [a X - ^ ^-^N + [a^ x • [at x + -G2{SO{7)) - -F^j , (54) 
[[at X ar X [h X h]r = I (^A^^ + - ^[^t x h]' ■ [at x ^C2(SO{7)) + , (55) 

[[at X at]*^ X [a X = -1= (^^N' - + ^[^^ x ■ [at x 5]^ + ^-^2(50(7)) + ^i^^) . (56) 

Using these relations, we can rewrite the Hamiltonian as 

^ = ^ V2F+lgF[[at X at]^ x [5 x 5]-^]" - 2V7p[a^ x 5]J 

F=0,2,4,6 



1 

m 
1 

2n 



eoN{N - 1) + ei(F2 - UN) + ^S+S^ + eg [at x at]^ . [5 x a 



e47V(/V + 6) + e5/^ + e6C'2(S'0(7)) H-er-F^ + e3[at X . [at X - pF,, (57) 
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where S+ = ^[d) x a^]°,5_ 



S\, and the couphng constants (z = 0, • • • ,7) are given by 



eo 

62 
64 
66 
67 



9 2 

21 10 



11 



54 



50 
1 



1 

YY56, 

18 7 
-.94 + 7736 



11- 



24 



1 



11- 
10 156 



^30 
1 

"7 
5 

'84 



11 



34 



^56, 65 — --50 



-32 



77 



■34 



262 



■96, 



90 



21 



32 



77 



34 



231 



36, 



32 



154 



34 



132 



36- 



(58) 



From Eq. (|57p we see that, if 63 — 0, dynamical symmetry alone can solve the problem, because the Hamiltonian 
would then be expressed in terms of the Casimir operators of the U{7), S0{7), 50(3), and 50(2) groups only. We 
shall henceforth discuss the problem in this situation. 

As in the case of spin-1 and spin-2 BECs, the quantum numbers of N,F, and arise from the U{7), 50(3), and 
50(2) groups. The quantum number of the 50(7) group is determined from the SU{1, 1) algebra: 



[5z, 5-t 



ih5^, [5^,5- 



-25„ 



(59) 



1 - 7 

where Sz = —N + —. As in the spin-2 case, the Casimir operator of the 50(7) group is written as 

02(50(7)) = N{N + 5)~4S+S- 
= 45^-^. 



(60) 



The eigenvalue of 0*2(50(7)) is 



02(50(7))|5,5, > 



45(5-1) 



21 



15,5, > 



= t{t + 5)\N,T > 



(61) 



This means that the quantum number of the 50(7) group is r. As a result, we have determined four quantum 
numbers; however, seven quantum numbers are needed to uniquely characterize the state. Therefore, there must be 
three missing labels, and the eigenstate is expressed as 

\N,T,r,q,s,F,Fz>, (62) 

where r, g, and s are missing labels [2^1 . The corresponding eigenvalue of the Hamiltonian is then given by 



1 
1 

20 



eoN{N - 1) + ^(iV _ t){N + t + 5) + ei{F{F + 1) - 127V} 
[eiN{N + 6) + e^N + eeT{T + 5) + ejFiF + 1)] - pF^. 



(63) 



The eigenspectrum is degenerate with respect to missing labels r, q, and s. 

VII. LOW-LYING EIGENSPECTRA AND EIGENSTATES 

In this section, we study the low-lying eigenspectra and eigenstates of spin-1 and spin-2 BECs which are valid for 
an arbitrary range of coupling constants. 



A. spin-1 case 



As can be seen from Eq. (PO)) . the ground-state phase is determined by the sign of ci only; it is ferromagnetic 
if ci < and antiferromagnetic if ci > 0. This result is consistent with the prediction of the mean-field theory. 
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The difference arises when the ground state is antiferromagnetic. In this case, while the mean-field study predicts 
F = 0, the exact result shows that F is or 1 according to whether the number of particles N is even or odd. The 
low-lying excitation spectra show this even-odd parity effect as illustrated in Fig. [TJ The F ^ 1 BEC of ^^Rb is 
ferromagnetic (ci < 0) and that of ^'^Na is antiferromagnetic. Their low-lying spectra are respectively illustrated in 

Fig. 1 (a) and (b) or (c). Since the ratio -i is of order 1, where g = — 'j^f~ ^'^^ as is the Bohr radius, a marked 

difference should arise in the low- lying energy- level spacing between the ferromagnetic and antiferromagnetic phases: 
the low-lying energy-level spacing of the ferromagnetic BEC is of the order of N, but that of the antiferromagnetic 
BEC is of the order 1 . This implies that there are numerous quasi-degenerate energy levels above the ground state of 
the antiferromagnetic BEC, and the ground state is therefore vulnerable to symmetry-breaking perturbations. 

B. spin-2 case 

The exact ground-state phases and low-lying eigenspectra of a spin-2 BEC are much richer than those of a spin- 
1 BEC since the number of the coupling constants and quantum numbers increases from 1 to 2 and from 3 to 5, 
respectively. The phase boundaries of the mean-field ground states of a spin-2 BEC have been shown as follows 

ferromagnetic : di < and + 4di < 0, (64) 
antiferromagnetic : > and ^5 4- idi > 0, (65) 
cycHc : di > and < 0. (66) 

In many-body theory, from the analysis of and (|47p . the ground-state phase boundaries are given as follows: 

4A^ + 2 

ferromagnetic : di < and ^5 + ^ ^ ^ di < {N ^ 2k, F ^ 2N, t ^ N, n^o ^ 0), (67) 
di < and d5 + ^^^T-fH^rif < {N = 2k + l, F = 2N, T = N, n^o = 0), 

(68) 

47V -I- 2 

antiferromagnetic : ds > and ds + ^ ^ ^ di > {N — 2k, F = 0, t = 0, 7130 = 0), (69) 

d5 - ^di > and ds + ^^j^ ~ j^^"^^^^ di > (TV = 2fc -t- 1, F = 2, r = 1, ngo = 0), 

(70) 

3 

di > and < ds < -di (iV = 2fc + 1, F = 0, t = 3, 7130 = 1), (71) 

N 

cyclic : di > and dg < (iV = 3fc, F = 0, t = N, n^o = — ), (72) 

o 

di > and ds + ^J^^^ di < (iV = 3fc + 1, F = 2, r - TV, ngo = (73) 

3 TV - 4 

di > and - -^j^^di < ds < (TV = 3fc + 1, F = 0, r = TV - 4, 11.30 = (74) 

3 TV — 2 

di > and ds + ^^^^ di < (TV = 3fc + 2, F = 2, r - A^, 1130 = )> (75) 

3 TV — 2 

di > and - 2N+l '^^ < 4 < (TV = 3fc + 2, F = 0, r = TV - 2, 7130 - )> (76) 

where fc S Z. Thus, the phase boundaries and their number change compared with the case of mean- field theory 
due to a finiteness of the number of particles. In fact, in the limit of TV — > 00, the phase boundaries between the 
ferromagnetic and antiferromagnetic phases reduce to those of mean-field theory. If TV = 6k, the exact ground-state 
phases have a one-to-one correspondence to those of mean- field theory. If we consider the cyclic phases with TV = 3fc + l 
or TV = 3A; -|- 2, the phases (|7i)) or (|75|) disappear and the cyclic phases correspond to the phases ([75)1 or ([75l) in this 
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limit. However, within the antiferromagnetic phases, the phase boundaries which appear with iV = 2fc + 1 do not 
disappear in this hmit; thus, the antiferromagnetic phases are always divided by ([70]) and ([7T|) . 

This leads to a marked difference with the spin-1 case because the low-lying excitation spectra of a spin-1 BEC 
are uniquely determined when the ground-state phases are specified. In Fig. [21 O and 21 we give diagrams of the 
low-lying eigenspectra up to the second excited states in the ferromagnetic phase with N — 2k, antiferromagnetic 
phase with N = 2k, and cyclic phases with N = 3k, respectively. The other particle-number cases can also be found 
in a similar way and a rich variety of cases appear in considering the higher excitation. In the ferromagnetic case, 
while the diagram (a) disappears in the limit N oo, the other diagrams do not disappear in this limit. Similarly, 
while in the antiferromagnetic phase all diagrams do not disappear, in the cyclic phase the only remaining diagram 

is (a) in this limit. Since — !- is of order of 1 ^ 10 for spin-2 species [1, [2^, [2^, the low-lying energy-level spacing 

of the antiferromagnetic and cyclic phases are of order of 1, whereas that of ferromagnetic phase is of order of N. 
Considering the huge degeneracies of the ground states found in the antiferromagnetic and cyclic phases [13], it is of 
interest to study the response of these ground states to symmetry-breaking perturbations. The investigation of this 
problem is underway. 



VIII. SUMMARY AND DISCUSSIONS 



In this paper, we have shown that dynamical symmetry completely determines the exact eigenspectra and eigenstates 
of spin-1 and spin-2 BECs in the single-mode approximation. In particular, a spin-2 BEC in this approximation 
corresponds to the C/(5) limit of the interacting boson model of atomic nuclei. We have also shown that dynamical 
symmetry alone cannot solve the eigenvalue problem in the spin-3 case, because the Casimir operator of the exceptional 
group G2 is proportional to that of the 5*0(7) group. We have, however, identified the class of the coupHng constants 
for which the exact eigenspectrum can be found by dynamical symmetry alone. 

Compared with the spin-1 case, the new term that appears in the spin-2 case is the term Vq which reflects the 
S0{5) symmetry in the spin-2 case. If we put di = p = in Eq. (I34p . we find that the symmetry of the Hamiltonian 
is not U{1) X S0{3) but U{1) x SO{5). We can generalize this fact to the spin-/ case because Vq is written as 



2/ + 1 



N{N + 2f-l)~C2iSO{2f + l)) . (77) 



That is to say, the Hamiltonian has the U{1) x 50(2/ + 1) symmetry within a specific class of coupling constants. 

The eigenspectrum of a spin-1 BEC in Eq. (I26p is characterized only by the total number of particles and the total 
angular momentum, which reflect the gauge invariance and the isotropy of space respectively. The eigenspectrum 
of a spin-2 BEC in Eq. (|34p involves an additional quantum number r, which does not reflect either the gauge or 
space-time symmetry but arises from the Casimir operator of the 50(5) group; r denotes the number of particles 
not in pairs of F = 0. The missing label n^o, which describes the number of spin-singlet trios, appears only in the 
eigenstate (|46|) and not in the eigenspectrum ((47l) of a spin-2 BEC. Within the specific class of the parameters with 
63 — 0, the eigenspectrum of a spin-3 BEC is completely characterized by the same set of quantum numbers as that 
of a spin-2 BEC, that is, N, F, Fz, and t, and it is in this class that the eigenvalue problem of a spin-3 BEC can be 
solved by the dynamical symmetry alone. 

We have also discussed the low-lying eigenspectra and eigenstates of spin-1 and spin-2 BECs in the absence of 
external magnetic field. For the spin-1 case, the phase boundary between the ground states with (p9)) and pO)) 
corresponds to that of mean-field theory. The exact low-lying excitation spectra are uniquely determined by the sign 
of Ci- On the other hand, for the spin-2 case, the phase boundaries between the ground states with (|46p and (|47p are 
different from those of the mean-field theory. These differences disappear in the limit N ^ 00 except for the phase 
boundary within the antiferromagnetic phase with iV = 2fc -I- 1. However, since all experiments of BECs in the cold 
atomic gases are done for a finite number of atoms, these differences seem to be real issues. In addition, the low-lying 
excitation spectra cannot be specified even if the ground states are specified, and the coupling-constant dependence 
of the low-lying excitation spectra is stronger than those of the ground states. This is because while the quantum 
number which means the number of particles not in pairs of _F = is equal to the quantum number F for a spin-1 
BEC, it is not so for the spin-2 case in which there arises an additional quantum number n^Q which connect these 
quantum numbers. We also show that for the case of spin-1 and spin-2 BECs, the low-lying energy-level spacings of 
antiferromagnetic and cyclic phases are by a factor of N smaller than those of the ferromagnetic phase. This means 
that for the antiferromagnetic and cyclic phases quasi-degenerate spectra emerge which may cause symmetry-breaking 
transitions to yet unexplored many-body states. 

Note added.- At the time of submission of this paper, we became aware of a paper [2^ by Van Isacker and Heinze 
who discuss the exact ground-state phase structure with arbitrary spin in the absence of external magnetic field. 
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APPENDIX: RECOUPLING FORMULAE 

The Wigner 6-j symbol is defined by [2] 



J Jl J2 I — \ ^ f-^\jl+j2+j3+mi+m2+m3+lll+ll2+fJ.3+l^l+1^2+l^3 

\ mi 7712 n^3 J ^—^ 

p,l,fJ.2p.3Ml,l^2M3 





j2 











Ml Ai3 



7l 72 73 \ 

where is the Wi gner 3-j symbol and is defined by 



^Ml Ms/ y/2j3 + 1 

The Wigner 6-j symbol is related to the Wigner 3-j symbol by [2] 



I mi 1712 1713 J \rn[ -msj f-^, V J3 7 



jl j2 jsl / jl j2 j'3 \ f jl j2 J3^ 
Jl J2 J3 J \™'l ™2 "^3/ V"^l "^2 -"T-3 



Using Eq. (|A.3p . one can derive the following recoupling formulae: 
For the spin-1 case, 



For the spin-2 case. 



I' 

For the spin-3 case. 



(A.l) 



- ^ ' <Jimj2M2|j3,-M3>. (A.2) 



(A.3) 



[at X a]' • [at X a]' ^ (2? + 1) ^ I J J [at x at]'' • [h x 5]'' + (A.4) 
at X ■ [at X a]' - {21 + 1)Y. [l 2 ^/ } ^^J'' ' ^l'' + ^^'^^ 



[at X ■ [at X 5]' = (2/ + 1) XI {3 3 \^ l^'' ^ ^^J'' ' ^ + (^-6) 
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FIG. 1: Low-lying eigenspectra of a spin-1 BEG. The angular momentum F of each state is shown to the right, (a) ci < 0. 
(b) ci > with even N . (c) Ci > with odd N . The low-lying energy-level spacing is of the order of N for the ferromagnetic 
BEG (ci < 0), while it is of the order of 1 for the antiferromagnetic BEG (ci > 0). /^p(^af) denotes the energy gap of the 
ferromagnetic (antiferromagnetic) ground state. 
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FIG. 2: Low-lying cigeuspcctra of a spin-2 BEG for the ferromagnetic phases with N = 2k and di < 0. The values of F, nao, 
and r of each state are shown to the left, middle, and right, respectively, (a) ds + f^^di < 0. (b) — ^^^j' rfi < ds < — ^j^~^ di. 



V^^ 2JV+1 

state. 



2JV+1 



di. (d) 



2N+1 



di < 



2JV+l"'i ^ \"/ 2N+1^ 2]V4 

< — ^^di. Af denotes the energy gap of the ferromagnetic ground 
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FIG. 3: Low-lying cigouspectra of a spin-2 BEG for the antiferromagnetic phases with N = 2k. The values of F, n^o, and r of 
each state are shown to the left, middle, and right, respectively, (a) di > and jdi < d[, < (b) di > and dr, > (c) 

di > and < ds < |di. (d) di < and ds > —^ff^di. Aaf denotes the energy gap of the antiferromagnetic ground state. 
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FIG. 4: Low- lying cigcnspectra of a spin-2 BEG for the cyclic phases with N = 3k and di > 0. The values of F, nso, and r 
of each state are shown to the left, middle, and right, respectively, (a) ds < — 2rJ+i di ■ (b) 



2JV+1 



.di. (d) 



-jw+idi < dr, < -2l^di- (c) 
- 2N-3 di < ds < 0. Ac denotes the energy gap of the cyclic ground state. 



